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LINEAR  THEORY  OF  MICROPOIAR  ELASTICITY 

A.  Cem&l  Eringen 
Purdue  University 


ABSTRACT: 

Equations  of  motion,  constitutive  equations  and  boundary 
conditions  are  presented  for  a  special  clasB  of  micro-elastic 
materials  called  Micropolar  Solids.  These  solids  respond  to 
micro-rotational  motions  and  spin  inertia  and  can  support 
couple  stress  and  distributed  body  couples.  The  couple  BtresB 
theory  is  shown  to  emanate  as  a  special  case  of  the  present 
theory  when  the  motion  is  constrained  so  that  micro-  and 
macro-rotations  coincide.  Several  energy  and  uniqueness 
theorems  are  given. 
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1 .  INTRODUCTION 

The  theory  of  micromorphic  materials  introduced  by  Eringen 
and  his  co-workers  [l]  to  [5]  deals  with  a  class  substance  which 
exhibits  certain  microscopic  effects  arising  from  the  local  struc¬ 
ture  and  micromotions  of  the  media.  These  materials  can  support 
stress  moments  and  body  moments  and  are  influenced  by  the  spin 
inertia.  The  general  theory  is  however  very  complicated  and  even 
in  the  case  of  constitutively  linear  elastic  solids  [2]  the  differ¬ 
ential  equations  are  not  easily  amenable  to  solution. 

To  simplify  the  theory  Eringen  and  Suhubi  [2]  introduced  a 
theory  of  couple  stress.  In  view  of  the  fact  that  in  the  abundant 
recent  literature  the  terminology  of  couple  stress  is  used  in  a 
different  context  we  thought  it  may  save  the  reader  from  confusion 
if  vs  name  our  theory  "Micropolar  elasticity".  In  contrast  to  the 
couple  stress  theory  in  micropolar  elasticity  all  components  of  the 
as8ymmetric  stress  tensor  are  determinate  and  the  motion  of  the 
media  is  fully  described  when  the  deformation  and  micro-ro*ation 
vectors  are  known.  The  concept  of  micro-rotation  and  the  corres¬ 
ponding  field  equations  are,  of  course,  totally  absent  in  the  couple 
stress  theory. 

Physically,  solids  that  are  composed  of  dipole  macromolecules 
may  be  adequately  represented  by  the  model  of  the  micropolar  elasticity. 
Fibrous  materials  and  some  granular  and  porous  solids  may  also  fall 
into  the  domain  of  this  theory. 

In  Art.  2  we  recapitulate  the  linear  theory  of  microelasticity 
given  in  [2]  briefly  in  slightly  different  form.  Art.  5  is  devoted 
to  the  derivations  of  the  basic  equations  and  boundary  conditions  of 
micropolar  elasticity.  In  Art.  4  it  is  shown  that  the  couple  stress 
theory  results  as  a  special  case  cf  the  micropolar  theory  when  the 
motion  is  constrained  so  that  micro-  and  macro -rotations  coincide. 


2 


Art.  5  is  devoted  to  a  discussion  of  energy.  Here  we  determine 
the  restrictions  that  must  be  imposed  on  the  material  constants 
in  order  that  the  internal  energy  be  non-negative.  Two  theorems 
on  uniqueness,  one  for  the  dynamical  solutions,  one  for  the 
static  case  are  given  in  Art.  6. 
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2.  BASIC  EQUATIONS 

The  theory  of  micro -elastic  solids  formulated  in  [1]  and  [2] 
are  based  on  a  set  of  lavs  of  motion  and  constitutive  equations 
some  which  are  new  to  the  mechanics  of  continue  and  others  are 
modifications  and  extensions  of  the  veil -known  principles.  The 
basic  lavs  are: 

Conservation  of  mass: 


t? ♦  (cV’k 


in  'V  (2.1) 


Conservation  of  micro-inertia,  [J]: 
^  i-km 


- —  +  i,  v-i  v  ,  -  i ,  v  =0 
km,r  r  rm  rk  kr  rm 

in  V 


(2.2) 


Balance  of  momentum: 


W  +  0(f<  -  V  ■  0 


in  Y  (2.3) 


Balance  of  first  stress  moments: 


*mi  8  mi  +  \im,k  +  p(/,  -  a  )  »  0 

£m  Jim 


in 


~V 


(2.4) 


Conservation  of  energy: 


oe=t  v  +(s  -t  )  v 
M  ki  l,k  v  ki  k \l 

+  \  ,  v  .  ,  +  q,  ,  +  ph  in  Y 

kirn  m/,k  %k 


(2.5) 


h 


Principle  of  entropy: 

pr=  pti  -  (^k  -  ^>0  in  1/  (2.6) 

Thie  inequality  is  being  axiomatized  to  be  valid  for  all 
independent  processes.  In  these  equations 


p  ■  mass  density 

v  ■  velocity  vector 

A 

i  *  i  «»  micro-inertia  moments 
km  mk 


k  / 


ki 


ki 

\im 

^im 

o, 

im 


*k 

h 

n 

e 


«  gyration  tensor 

■  stress  tensor 

*  body  force  per  unit  mass 

*  8  *  micro-stress  average 

L  K 

■  the  first  stress  moments 

■  the  first  body  moments  per  unit  mass 

■  inertial  spin 

=  internal  energy  density  per  unit  mass 
=  heat  vector  directed  outward  of  the  body 

»  heat  source  per  unit  mass 

*  entropy  per  unit  mass 

■  temperature 


Throughout  the  present  paper  a  rectangular  coordinate  system 
(V  x2,  x5)  is  employed,  Pig.  1.  An  index  followed  by  a  comma 
represents  partial  differentiation  with  respect  to  space  variables 
and  a  superposed  dot  indicates  material  differentiation,  i.e. 
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k ,/ 


(2.7) 


Here  and  throughout  this  paper  repeated  indices  denote  summation 
over  the  range  (1,2,3).  The  foregoing  equations  are  expressed  in 
Eulerian  representation.  However,  in  the  linear  theory  of  solids 
the  difference  between  the  Eulerian  and  material  representations 
disappears  since  in  calculating  the  time  rates  the  convective 
terms  are  ignored.  Expressions  (2.1)  to  (2.6)  are  valid  at  all 
parts  of  the  body  B  having  volume  '‘V  and  surface  ,  except 

at  finite  number  of  discontinuity  surfaces,  lines  and  points.  At 
the  surface  S ^  of  the  body  we  have  the  boundary  conditions 


nk  3  ^n)! 


\im  k  ^(g)/m 


on  y  (2-8) 

on  y  (2.9) 


where  n  is  the  exterior  normal  to  and  t.  ,  and  A.  .  are 

J  (a)  "(a) 

respectively  the  surface  tractions  and  surface  moments  acting  on 
^  .  For  the  spin  inertia  we  have  the  kinematical  relation 


o,.  =  i  .(v  +  v  v  ) 

kz  ml  mk  nk  mn 


(2.10) 


We  note  that  while  equations  (2.1)  and  (2.3)  are  well-known 

from  the  classical  theory,  six  equations  (2.2)  are  entirely  new 

and  nine  equations  (2.4)  and  the  energy  equation  (2.5)  are  modified 

and  extended  forms  of  the  corresponding  limited  axioms.  In  fact 

the  skew-symmetric  part  of  (2.4)  is  certainly  new  and,  of  course, 

the  kinematical  variables  v  „  and  i,  ,  do  not  enter  into  the  classical 

kz  kz 


continuum  mechanics . 
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If  we  exclude  the  phenomena  of  heat  conduction,  in  the 
present  theory,  the  determination  of  motion  requires  the  deter 
mination  of  the  following  nineteen  unknowns 


p(*>t)  ,  ^(^t)  ,  vk(2,t)  ,  vki(2,t)  (2.11) 

as  against  the  four  unknowns  p  and  v  of  the  classical  theory. 

k 

The  basic  equations  (2.1)  to  (2.5)  are  valid  for  all  types  of 
media  independent  of  their  constitutions.  The  character  of  the 
media  are  reflected  through  the  constitutive  equations.  For  micro¬ 
elastic  solid  in  [2]  we  gave  sets  of  constitutive  equations  valid 
for  nonlinear  deformations .  Here  produce  only  the  linear  theory 
relevant  to  the  present  work.  These  are 


t  =  [(A  +  r)e  +  ri  e  )  6  +  2(u  +  o  )e 

k  l  1  rr  1  rr  k  t  11 


k  l 


(2.12) 


+  k  c  +  v  e 
1  fk  1  kf 


\t  ■  ((\  *  2T>err  +  <="1  '  ,)err>  6k,  +  2(|il  +  S°lK l 

+  (V1  *  K1  ■  °1>  (tw  *  V 


(2.15) 


\/m  “  Tl(\r,r  6mi  +  9rr,i  V  +  T2(q>rk,r  &mi 


+  ^rr,ro  \l ^  +  ^rr,k  &mf  +  T4  ^/r,r  &mk 


+  t(cd  6  +©  S  )  +  cp  6 

5'  rl,r  mk  mr,r  kt  6  rm,r  ki 


(2.14) 


+  T  (p  +  T  (CP  +  CD  )  +  T  CD 

7  *m,k  8'  mk,l  Kl,in  9  *k,m 


+  T  tp  +  T  <? 

10  mi,k  11  km,  .2 
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where  X^,  x,  v  and  to  are  material 

constants  and  e,  .  and  e  ,  are  respectively  the  strain  tensor  and 
kl  k  i 

the  microstrain  tensor  of  the  linear  theory  defined  by 


ew  si(V  +  V) 


(2.15) 


"ki  *  \,t  + 


(2.16) 


It  can  be  shown  that  [2]  the  arc  length  and  the  single  changes 
can  be  calculated  once  these  quantities  and  ffi  are  known'*’. 

y  IQ 

For  the  linear  theory  we  also  have 


v  ■  © 
ki  vki 


(2.17) 


A  thermodynamic  discussion  can  be  made  showing  that  for  a  non¬ 
heat  conducting  media 


q  ■  0 


*1  ■ 


de 


(2.18) 


-l 


where  \|r  =  e  -  is  the  free  energy  which  depends  on  p  and  6 
only.  Alternatively 


v  and  9  used  here  correspond  to  v  and  qp  of  [1]  and  [2] 
ki  ki  ia  ik 

and  v  .  and  (p  .  of  [j]  smd  [5J. 
ki  ki 

2 

For  these  and  other  interesting  results  we  refer  the  reader  to 


[1],  [5i  and  [5] . 
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q  =  0  ,  e  =  r-  (2.19) 

-  OT)  P 

where  e  is  considered  to  be  a  function  of  t)  and  p  \ 

Upon  substitution  of  (2.15)  and  (2.16)  into  (2.12)  and  (2.13) 
and  the  result  and  ( 2 . 14 )  into  (2.3)  and  (2.4)  together  with  (2.1), 
(2.2),  (2.10),  (2.17)  and  v  =  £  we  obtain  just  the  right  number 
of  partial  differential  equations  to  determine  the  nineteen  unknowns 

p(*>t)  ,  ^(*>4)  ,  uk(*,t)  and  <Pk/U>t). 

Under  appropriate  boundary  and  initial  conditions  these  partial 
differential  equations  should  be  adequate  to  predict  the  behavior 
of  micro-elastic  solids . 

The  theory  outlined  above  is  too  complicated  for  application. 
The  theory  of  raicropolar  elasticity  presented  below,  however,  not 
only  brings  some  elegance  into  the  picture  but  also  is  simple 
enough  for  mathematical  treatments . 
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3.  THEORY  OF  MICROPOLAR  ELASTICITY 


Definition.  A  body  vill  be  called  micropolar  elastic  solid 
if  for  all  motions. 


kmi  ’ 


9,.,  -  9. 


(5.1) 


The  motion  of  micropolar  solids  appears  to  have  important 
implications.  Such  materials  are  physically  realistic.  They 
exhibit  micro-rotational  effects  and  can  support  body  and  surface 
couples.  Materials  consisting  of  dumbbell  macromolecules  belong 
to  this  class.  Some  fibrous  and  granular  media  are  also  approxi¬ 
mated  by  this  model. 

By  use  of  the  skew  symmetric  properties  of  A  and  (£  the  basic 
equations  of  such  continua  can  be  simplified  a  great  deal.  We 
first  define  couple  stress  tensor  m^,  micro-rotation  vector  q>^, 
micro -angular  velocity  Vj  micro-rotation  inertia  a ^  and  body 
couple  1 b£ 

m  5  -  €  ,  \  . 

kr  rim  klm 


©  =  —  €  © 
r  2  rkl 


v  b  —  e  v 
r  ”  2  rkl  kl 


(5.2) 


a  =  -  e  c 
r  “  rkl  kl 


1  =  -  6  1 
r  rkl  kl 


The  couple  stress  tensor  m  has  the  same  sign  convention  as  the 

kr 

stress  tensor  t  ,  Fig.  2. 
kr 

Upon  using  (3.2)  in  (2.4)  and  (2.5)  these  equations  can  be 
expressed  respectively  as  [5). 
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m  +  €  t  +  p(/  -  a  )  m  0 

rk,r  kir  It  v  k  k' 


(3.5) 


PE  '  \l(vi,k  ’  eMr  V  +  "k*  V 

+  \,k  +  0h 


(3.4) 


The  boundary  conditions  (2.9)  can  be  written  as 


m  n  =  m.  . 
rk  r  (u.)k 


on 


S' 


(3-5) 


where  m^  ^  is  the  surface  couple  acting  on  ^  . 

Similarly  the  constitutive  equations  for  t  and  \  can  be 

reduced  to  simpler  forms.  In  simplifying  the  expressions  of 

\  we  make  use  of  the  skew-symmetry  property  of  \  and  $  as 

expressed  by  (3-1)  and  employ  the  expressions  (3.2) .  This 

leads  to  the  following  restrictions  on  the  coefficients  t  . 

a 


T1  "  T2 


\  bT6 


T  =  T 

9  11 


(3.6) 


Using  (3.2)1  and  (3*2)^  and  introducing  a  set  of  new  elastic 
coefficients  by 


a  -  2(t9  -  Tg)  , 


6  £-2(t4  -  T^) 


y  £  2(-t|  +  T  -  T  +  TQ  -  T  +  T,rt) 

“  4  5  7  3  9  10 


(3.7) 


A  £  ^  +  T  +  rj  ,  \1  b  |j  +  al  +  v1 


*  S  *1  “  V1 


(3.8) 
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*kl  ■ \t  <5-13> 

Carrying  this  into  (2.2)  and  recalling  that  v  =  -  v 

&.Z  Z& 

we  obtain 
Di 

—  »  0  or  i  b  constant  =  j/2  on  material  lines  (3.1'0 

Dt 

With  this  then  the  expression  of  inertial  rotation  becomes 


(3.15) 
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The  field  equations  of  linear  micropolar  elasticity  are 
obtained  by  substituting  (3*9 )>  (3*10)  into  (2.3)  and  (3>3)> 
Hence1 


(X  -  +  <P  ♦  K\>tl  +  +  P(ft  -  \)  -  0  (3.16) 


} 


(a  +  |3,*Mk  + 


yep  +  Ke  u 

kfm  m,i 


2x\  +  =  0  (3.17) 


In  the  linear  theory  the  density  p  is  to  be  treated  as  constant 
so  that  no  need  arises  for  the  use  of  the  continuity  equation  (2-1). 
Generally  the  integral  of  this  equation  for  the  linear  theory  is 
used  for  the  purpose  of  calculating  small  changes  in  density  after 
the  displacement  field  is  known.  This  is  given  by 


p°/p  E 1  +  \,k 


(3-18) 


The  energy  equation  (3*4)  with  the  use  of  (3*9)  and  (3*10) 
becomes 


*  ■  Xekkdi/ +  (2p  +  "’v1* +  2*(rk  -  V(,"k  'V 


-  “V  \,t +  p  W +  k 


(3.19 


\ 

) 


where  w  is  the  vorticity  vector  defined  by 

A 

1 

U  -  —  €  V 

k  2  kim  m,l 

Equation  (3-19)  can  be  integrated  with  respect  to  time  to  give 


(3-20) 


’These  equations  can  be  made  to  coincide  respectively  with  eqs.  (6.13) 
and  (S.l4)  of  [2]  by  redefining  material  constants  and  using  (3.2)2- 


15 


<*  -  po£o  -  5  !x  V  eu  *  (2p  +  K)  \i  eik 


♦  2*<  V  \)(rk  -  V  +  “V  %t 

♦  Mk,/  v +  >v  v 


(3.21) 


where  pq  and  are  the  mass  density  and  internal  energy  in  the 

undeformed  body.  Without  loss  in  generality  we  may  take  e  =  0. 

o 

In  obtaining  (3-21)  in  the  spirit  of  the  linear  theory  we  assumed 
that  g,  =  An  alternative  form  to  (3 .21)  is 


*  *  i  +  'fcta  '  V’  **U*t,*i 


(3.22) 
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4.  COUPLE  STRESS  THEORY 


We  now  examine  a  special  case  of  the  above  theory  namely 
the  motion  with  constraints 


<P,. 


(4.1) 


In  this  case  the  stress  constitutive  equations  (3*3)  reduce  to 


Xe 

rr 


+  (2m  +  K)eu 


which  are  identical  to  those  of  the  classical  theory  of  elasticity 

K 

provided  we  replace  p  +  g-  by  p.  The  field  equations  (3.16) 
now  become 


(x  + 11  + i>  V k +  ^  *  f  \ ,it +  p(fk  •  \> 


which  are  the  celebrated  equations  of  Navier  with  the  replacement 
of  p  for  p  +  k/ 2.  From  equations  (3*17)  u  drops  out.  Since  q> 
is  no  longer  an  independent  variable  (3.17)  fully  determines  the 
body  couple  £  admissible  for  this  class  of  motions.  With  this 
viewpoint  then  the  present  theory  goes  into  the  classical  theory 
of  linear  elasticity. 

There  is,  however,  another  way  of  interpretation.  Suppose 
that  (4.2)  only  determines  the  symmetric  part  of  the  stress  tensor 
i.e.,  write 


rr 


+  (2p  +  Kje^ 


(4.2) 


15 

where  a  (  )  enclosing  the  paranthesis  indicates  the  symmetric 
pare  of  a  tensor  and  similarly  a  [  ]  will  be  used  to  indicate 
the  skev-symmetric  part  of  a  tensor,  i.e. 


a(ki)=  2  (aki  +  V 


[ki] 


*  2  '  aik^ 


We  now  solve  for  the  antisymmetric  part  t^^  of  the  stress 

tensor  from  the  original  equation  (3*3)  of  the  balance  of  moment 

of  momentum.  By  multiplying  (3*3)  by  €  we  first  put  it  into 

kmn 

the  form 


tr  +  -  e  „  m  -  p(ir  -  or,  .,)  **  0 
[ki]  2  rki  nr,n  [ki]  [ki]y 


(4.3) 


Upon  substituting  (4.1)  into  (3*10)  and  (3*15)  with  (3*12)  used 


for  r  we  get 

it 
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m  -  -  e  .  u 
ki  2  imn  n,mk 


(4.4) 


4tki]  ■  J  (4-5) 

Substituting  these  into  (4.3)  we  find 

‘[ki]  '  2  (\,inn  •  “/.km,1  +  P<*[ki)  ’  J  U[k,i]>  <k-6) 

Using  (4.2)  and  (4.6)  we  calculate 

\i,k  *  ^(ki),k  T  ^[ki^k 


and  substitute  it  into  the  equations  of  motion  (2.3) •  Hence 


l6 


0,+  ^+  i+  Iv2).  +  (n  +  £  -  l  V^)u 

x  2  4  '  k,ik.  N  2  4  '  i,kk 

+  7is\-.  J  +pf1B  p(1-  i  y2)  “i 


2  ki,k  /k,k 


(4.7) 


Pj  •• 

+  —*■11 

2  k,/k 


where  V  is  the  laplacian  operator  in  rectangular  coordinates,  i.e. 


=  u 

"  i,kk 


Equations  (4.7)  are  the  field  equations  of  a  theory  jcnovn  as  the 
couple  stress  theory.  The  present  theory,  therefore,  for  a  con¬ 
strained  motion  gives  the  couple  stress  theory. 
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5-  NON -NEGATIVE  INTERNAL  ENERGY 

In  this  section  we  investigate  the  conditions  that  the  elastic 
moduli  must  be  subjected  to  in  order  that  the  internal  energy 
density  be  non-negative.  These  conditions  have  important  impli¬ 
cations  in  regard  to  stability  problems,  wave  propagations  and 
uniqueness  theorems . 

Tneorem  1.  (non -negative  internal  energy) .  The  necessary  and 

sufficient  conditions  for  the  internal  energy  to  be  non -negative 
are 


0  <  3>\  +  2p  +  k 
0  <  3a  +  27 

The  sufficiency  of  (5*1)  is  proven  by  merely  observing  that  when 
these  inequalities  hold  each  one  of  the  following  energies  con¬ 
stituting  (3.21;  are  non-negative 


0  <  p  ,  0  <  K 

-7<P<7  ,  0  <  7 


(51) 


WE  2  i  [Xekke«  +  (2M  +  *>  \te, k1 


(5.2) 


2  K(rk  •  V  (rk  •  V 


2  i  [aw  +  + 


(5.5) 

(5-4) 


resulting  in 

p€  -  p(€e  +  eR  +  eM)  >  0  (5.5) 

The  fact  that  pc  is  non-negative  under  the  conditions  (5.1)  , 

E  1 

(J.l)^  are  well-known  from  the  classical  theory  of  elasticity. 
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In  fact  these  conditions  are  also  sufficient  to  make  pe  non- 

L 

negative.  It  is  simple  to  observe  that  for  r  *  cp  ,  pe  >  0 

K  K  ia 

if  and  only  if  k  >  0.  To  see  the  same  for  pe  we  write  this 

M 

expression  in  the  forms. 


p€u  "  7(3  a  +  M  7)(<P,  n  +  ^  0  +  9,  7)2  +  f>  (($,  +  %  J2 

M  3  1,1  2,2  3>3  1>2  2,1 

+  (9,  ,  +  9,  J2  +  ,  +  9,  J2  J  (r  +  P)  (9  2  +  9?  .  (5-6) 

1,3  3,1  2,3  -  3,2  1,2  2,1 

.2  -2  -2  _2  , 

+  cp  +  cp  +  cp  ,+  cp  1 

yl,3  t,l  2,3  3,2 


where  upper  signs  go  together  and  lower  together  and 


9 


r,r 


(5.7) 


From  (4.6)  we  see  that  in  either  case 


3a  +  6  +  y 

>  0 

3  >  0 

7-  0  >  0 

(5.8) 

3a  +  0  +  7 

>  0 

3  <  o 

7  +  j3  >  0 

(5-9) 

we  have  pe  >  0  so  that  (3.1)  are  sufficient  to  make  pe  >  0. 

M  —  M  — 

Conditions  (5-1)  are  also  necessary  for  the  nor.-negative- 

ness  of  pe .  To  see  this  we  first  note  that  e  .  r  -  cp  and  cp 

k  £  k  k  k,i! 

can  be  varied  independently  of  each  other.  Since  the  above  three 
energies  (5.2)  to  (5.4)  are  uncoupled  with  respect  to  these  vari¬ 
ables  each  one  of  these  energies  must  be  non-negative  independent 

of  each  ether.  Certainly  k  >  0  is  also  necessary  for  pe  to  be 

/  1  R 

non-negativs  provided,  of  course,  rj  cp  .  From  classical  elas- 

K  K 

ticity  theory  we  also  know  that  (6.1)^  and  (6.1)^  are  necessary 

and  suff  <'  ent  conditions  for  pe  >0.  Hence  we  would  like  to 

E  - 

"'"For  r  =9  (the  couple  stress  theory),  k  need  not  be  restricted  by 
k  k 

(5»l)y  However,  in  this  case  by  writing  2p  for  2p  +  K,  k  disappears 
from  all  equations. 
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find  out  if  the  last  three  conditions  in  (6.1)  are  necessary  for 

pe  to  be  non-negative.  To  this  end  vrite  (4.4)  as  a  quadratic 
M 

form  in  a  nine  dimensional  space,  i.e. 


0€  ■  a  v  y  .  a  =  a 

ij  yj  '  ij  Ji 


where 


h-  V 

1 

^2  ~  ^2,2 

’  y3 

2  ”3,3 

\  -  *1,2 

> 

in 

w8 

,  y6 

2  "V.,3 

y7  *  ^3,2 

) 

y8  ■-  ”5,l 

'  y9 

SS  CD 

*  1,3 

a  «a 

*  a  ,  « 

a  +  p  +  y 

a,  ■ 

a,,  -  a 

11  22 

33 

12 

13 

a45  “  a67 

s  m 

8  9 

P 

a, ,  ■  a 

b  a, .  * 

a  ■  a_„  ■ 

a  *  7 

44  55 

66 

77  88 

99 

all  other 

a  ■  0 

ij 

The  characteristic  values  a  of  a  are  obtained  by  solving  the 

1  IJ 

equation 

det  (‘ij  •  a8ij)  '  0 
The  nine  roots  a^  of  this  equation  are 


ai  "  a2  "  a3  *  7  "  P  ' 


a5  -  a6  “  a7  -  a8  "  7  +  P 


a^  =  yx  +  p  +  7 
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In  order  a  y  y  >  0  to  be  satisfied  to  all  y  it  is 
U  i  j  ~  i 

necessary  (and  sufficient)  that 


7  -  P  >  0  ,  7  +  p  >  0 

3a  +  p  +  7  >  0 

which  is  the  proof  of  the  theorem. 
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6.  THE  UNIQUENESS 


The  uniqueness  theorem  is  a  statement  that  the  field 
equations  (3.16),  (3*17)  or  in  curvilinear  coordinates 


(X  +  2|i  +  k)  2  2*  )i  -  (11  +  k)  2*2*H+k2*,P 
+  p(£  -  *i)  -  0 


(a  +  p  +  7)22-H-72*2*®+K2*li-2K:<p 


(6.1) 


+  p(£  -  J$)  -  0 


(6.2) 


possess  unique  solutions  under  certain  boundary  and  initial 
conditions  for  certain  values  of  elastic  moduli! . 

As  the  initial  conditions  ve  consider 

3lU,0)  -  l^U)  in  *Y  (6.3) 

JA(&0)  -  ^(jc) 

9(2^0)  -  <P0U)  in  *V  (6.4) 

j(x,0)  -  ^(x) 

and  are  prescribed  in  '~V  .  We  leave  the 
boundary  conditions  presently  unspecified.  Ultimately  we  intend 
to  show  that  variety  of  conditions  are  possible.  For  example  11 
and  <p  may  be  specified  on  the  boundary  surface  of  the  body. 

It  is  equally  permissible  to  specify  the  tractions  and  couples 
on  provided  that  they  are  in  equilibrium.  Many  mixtures  of 


where  V  V  ^ 


both  are  also  possible.  In  fact,  all  admissible  boundary  conditions 
allowing  unique  solutions  must  satisfy 
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*(a)k  \  *  5(a)k  \  "  0 


on 


y7 ,  t  >  o  (6.5) 

where  £,<?,£,  .  and  m.  .  are  respectively  the  differences  of 

-  (a)  (a) 

li,  <£,  and  from  their  respective  values  on  j  . 

Theorem  2.  (Uniqueness) .  Let  the  conditions  (5.1),  r  f  <p  and 

1  ^ 

j  >  0  be  satisfied  in  a  bounded,  regular  domain'  '"[/  of  space 
with  boundary  .  Then  there  exists  at  most  one  a(a>'t)  and 
one  j(;>c,t)  both  twice  continuously  differentiable  for  £  in  '1/* 

+  5^  and  0  <  t  <  w,  which  satisfy  (6.1)  and  (6.2),  the  initial 
conditions  (6.3),  (6.4)  and  a  set  of  boundary  conditions  compatible 
with  (6.5)  . 

To  prove  this  theorem  suppose  that  contrary  is  valid  and  two 
solutions 


(a)  (a)  ,  , 

a  ,  Z  *  (a  -  1,  2) 

exists  satisfying  (6.1)  to  (6.5).  Let 


(6.6) 


(1)  (2) 

a  £  a  -  a  > 


m  _(D  _(2) 

cp  1=  cp 


Then  clearly  u  and  q>  satisfy  (6.1)  and  (6.2)  with  iy> 
the  initial  conditions 


(6.7) 

i  *  0  and 


a(x>o)  =  ate;0)  =  cp(*,o)  =  <pU,o)  =  o 

•*/ 


(6.8) 


Let  t^  ,  ra^  ,  t^njk  and  m^  ^  te  respectively  the  stress,  couple 
stress,  surface  traction  and  surface  couple  corresponding  to  the 
solution  (6.7),  i.e. 


The  terminology  "regular  domain  of  space"  is  used  to  denote  domains 
for  which  the  Green-Gauss  theorem  is  valid. 
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m  (1) 
k t  ~  k£ 


*  •  • 


(6.9) 


We  compute  the  rate  of  the  total  internal  energy 


5 


—  /  pe  dv  = 


r 


dt , 


pe  dv 


"V 


where  we  employed  the  continuity  equation  -- —  =  0  .  Using 

(3-4)  with  q  =  0,  h  =  0  ,  this  may  be  put  into  the  form 


^  ■  W \  '  Wk/r 

t(\<V,k  '  'Ye,kV<JV 


Applying  Green-Gauss  theorem  to  the  first  and  fourth  term  and 
replacing  t  and  ra  by  their  respective  values  (eqs.  2.3  and 

KX  y  ic  KL  y  K 

(3.3)  viOi  t  =  1  =  0)  i.e. 


k£,k 


-  t  +  p j  v 

£kr  kr  £ 


we  obtain 


y 


(t  V  +  m,  v  )n  da 
'  k  £  £  ki  r  k 


t  +  X 

where  £  is  the  exterior  normal  to  and 

X  -l/p(vk\  + JVk,dT 


(6.10) 


(6.11) 


Eyi 


is  the  total  kinetic  energy .  If  the  boundary  conditions  are  such 
as  to  make  the  right  hand  side  of  (6.10)  vanish  we  get 


£  (t)  +X(t>  ■  c 
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where  C  is  a  constant, 
we  have 


Since  both 


t 


are  non -negative 


$  (0)  +X<0)  "  0 

Consequently 

tr  (t)  +  X  (t)  -  0  ,  0  <  t  <  oo  (6.12) 

With  j  >  0  we  have  ^  non-negative.  By  virtue  of  (4.1)  ^  is 

non-negative.  Therefore  we  must  have 

£  (t)  =  o  ,  3<f(t)  =  o  o  <  t  <  oo 


From  the  second  of  these  it  follows  that  since  (jc,0)  ■  <^£,0)  =  0 
it  follows  that 

u(x,t)  =  0  ,  <g(£,t)  in'V,  0  <  t  <  oo 

This  means  that  and  s  Hence  the  proof  of 

the  theorem. 

Next  we  investigate  the  static  case.  For  this  we  give  the 
following  theorem  which  is  the  counterpart  of  the  Clapeyron's 
theorem  in  the  classical  elasticity  theory. 


Theorem  5 .  If  a  body  is  in  equilibrium  under  a  given  system  of 
body  loads  £,  £  and  surface  loads  ^  and  jn^,  then  the  internal 
energy  of  deformation  is  equal  to  one  -half  of  the  work  of  external 
loads  that  produce  displacements  and  rotations  ^  and  qp  from  the 
unstressed  state .  Mathematically 
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/ 0(fk\  *  Vk}  dv  + 

V 

■  2  pc  dv 

Ar 


+  "(a)kV  da 


(6.13) 


The  proof  of  this  theorem  follows  by  forming  the  first  integral 
on  the  left  and  using  Green-Gauss  theorem.  For 

T  !  -/  c(fkV  W  dv  "  /“i k,A  +  (“/m 

v  cy 

+  Wlr  V  dv  ‘/( Vk  +  "£k  \\l  dv 

■/  (tkA,i  +  Vi,l  +  £klrt/A)dV 

*y 

By  use  of  Green-Gauss  theorem  we  convert  the  first  integral  on 
the  right  to  a  surface  integral.  Upon  using  (3-22)  in  the  second 
we  get 


v  -fy  (t(a)“k  +  “(aJk  V  da  20E  dv 

which  is  the  proof  of  the  theorem. 

The  above  theorem  can  be  used  to  prove: 

Theorem  4.  (Uniqueness  of  static  solutions) .  Let  the  conditions 

(5*1)  with  u  «  cp  ■  0  and  <p,  jt  r,  be  satisfied  in  a  bounded,  regular 

domain  *Y  of  space  with  boundary  y’  .  Then  there  exists  at 

most  one  stress  field  t.  and  one  couple  stress  field  m.  .  arising 
r  " "  lc if  1  — — .  "  k# 

from  twice  continuously  differentiable  displacement  field  ]a(2c)  and 
micro-rotation  field  <p(£)  satisfying 


4 
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t.  .  ■  t  ,  m.  .  •  m 

(&)k  ok  (ii)k  ok 


00  X 


(6.14) 


u  »  u  .  cp  *  cp 

k  ok  '  \  yok 


on 


y-y 


(6.15) 


Proof:  Let  (jc)  and  9^  (^c)  ,  (a  =  1,2)  be  two  distinct 
solutions  of  (6.1)  and  (6. 2) (with  u  e  9  e  0)  satisfying  (6.14) 
and  (6.15).  In  view  of  the  linearity  of  (6.1)  and  (6.2) 


(1)  (2) 
li  —  3i  "34 


m  m(l)  (2) 

>  Cp  E  9  -  9 


are  also  the  solution  of  these  equations  with  vanishing  £  and  £. 
Consequently  by  the  theorem  3  above  we  have 


l  (t(a)4\  +  B(a)4  V  06  dv 

But  on  y*  we  have  t.  .  ■  m.  .  *  0  since  both  tf  \  and 

L  (a)  (a)(!)  (2)  ia)  Aa) 

have  the  common  value  and  m|^j  and  m^  the  same  value  jn^. 

Similarly  on  ^  ^  ■  ?  »  0  .  Thus  (6.l6)  reduces  to 


(6.16) 


f  <* 


dv  =  0 


But  because  of  the  conditions  (5.1)  this  integral  can  vanish  only 
if 


%iu0  ’  \  '  rk  ■ 0 


(6.17) 


According  to  constitutive  equations  (3-10)  and  (5.11)  then  t. 

m  =0.  Hence  the  proof  of  the  theorem. 
k£ 


k  l 


The  integral  of  (6.17)^is 


\t  x< +  bk 


(6.:8) 
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where  fl  =  -  n  and  b  are  constants.  Substituting  (6.l£) 
k£  ik  k 

into  (6.17Lto  calculate  r  and  using  (J.12)  we  get 
L  k 


where 


CD  =  0. 

yk  k 


_  1 

\  2  Cyjn 


is  a  rigid  rotation.  We  therefore  see  that  the  displacement  and 
micro-rotation  fields  are  determined  within  a  rigid  rotation  and 
translation.  Wheny7^  ^  ^  0  we  have  the  mixed  boundary  conditions, 
Since  in  this  case  u  =  9  >  0  on  f/  -  ^  vc  see  that  «  b  *  0 

and  the  fields  ^  and  cp  are  uniquely  determined. 

In  the  case  we  have  the  traction  boundary  problem . 

For  this  case  the  fields  and  cp  are  determined  within  a  rigid 
rotation  and  translation.  When  to  and  ihq  are  in  equilibrium,  i.e. 

(f  t  .  dft  =  0  ,  <f  (x  *  ^  +  Hq)  .  da  =  0  (6.19) 

y  y 

then  u  and  cp  are  again  uniquely  determined.^  "Mixed-Mj.xed"  type 
of  boundary  condition  exists  for  which  the  solution  of  (6.1)  and 
(6.2)  is  unique.  All  such  class  of  solutions  must  satisfy 


t,  \,  u  +  m  cp,  =0 

(fi)k  k  (u)k  k 

where  ®(ji)  *  ^  an(*  ^  are  <***’*’erence  °f  ^n)  *  ®(]^),  ^ 

and  9  from  their  respective  values  on  y  . 
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